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Abstract 

We lay down the foundations of particle dynamics in mechanical theories that sat- 
isfy the relativity principle and whose kinematics can be formulated employing reference 
frames of the type usually adopted in special relativity. Such mechanics allow for the 
presence of anisotropy, both conventional (due to non-standard synchronisation proto- 
cols) and real (leading to detectable chronogeometrical effects, independent of the choice 
of synchronisation). We give a general method for finding the fundamental dynamical 
quantities (Lagrangian, energy and momentum) and write their explicit expression in all 
the kinematics compatible with the basic requirements. We also write the corresponding 
dispersion relations and outline a formulation of these theories in terms of a pseudo- 
Finslerian spacetime geometry. Although the treatment is restricted to the case of one 
spatial dimension, an extension to three dimensions is almost straightforward. 
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1 Introduction 

The goal of the present paper is to lay down the foundations of dynamics for mechanical theories 
satisfying the principle of relativity. 

We have shown elsewhere [TJ [2 J that, starting from the composition law for velocities and 
using the principle of relativity, the usual definition of kinetic energy for a particle (as a scalar 
quantity whose change equals the work done on the particle), and assuming the existence 
of elastic collisions between asymptotically free particles, one can construct such quantities as 
momentum, kinetic energy, the Lagrangian and the Hamiltonian for a free particle in an inertial 
frame — that is, all the basic ingredients one needs in order to build up dynamics. Although 
the treatment in Ref. pQ was restricted only to the pedagogically relevant cases of Newtonian 
and Einstein's mechanics, the basic technique is general and can be applied to produce a wider 
class of theories — essentially, all those in which the principle of relativity holds. 

How general are such theories? Galilei's and Einstein's kinematics require homogeneity of 
space and time, isotropy of space, the relativity principle, and a pre-causality condition (see 
Sec. 12.11 below for more details) o At first then, one might expect that one could generalise 
mechanics a good deal, by relaxing one or more of these hypotheses. However, on closer 
inspection it turns out that there is not much room left. Indeed, homogeneity is crucial in 
order to set up the very notion of a reference frame of the type commonly considered in special 
relativity. Moreover, as we wrote at the beginning, we are interested in theories that do satisfy 
the relativity principle. And, finally, pre-causality is fundamental if we want to be able to do 
physics at all. Hence, the only hypothesis we can drop is the one of isotropy. By this, we do not 
mean geometric isotropy (we still assume that space be Euclidean in every inertial frame) but 
the different notion of mechanical isotropy — that all oriented directions in space are equivalent 
in kinematics and free particle dynamics. 

The most general kinematics satisfying all the other postulates except mechanical isotropy 
were studied by Lalan in 1937 [3J, and will be reviewed and commented in Sec. |2] in order 
to provide a self-contained presentation. In these kinematics, anisotropy manifests itself in 
two ways. First of all, there could be different invariant speeds along different directions. 
Thus, if we assume that light in vacuum travels at the invariant speed, as usual0 it follows 
that in these theories it does not propagate isotropically. Furthermore, the factor that relates 
the measurements of time intervals by clocks in relative motion, and the analogous factor for 
measurements of lengths, depend on the direction of the relative velocity. Interestingly, these 
two manifestations of mechanical anisotropy are independent of each other — a circumstance 
which should not be entirely surprising, because the first type of effect is an artifact due to 
a non-standard convention for synchronisation of distant clocks, whereas the second one is 
objective and cannot be gauged away. In Sec. Owe describe a general method for constructing 
the fundamental dynamical quantities corresponding to a given relativistic kinematics, which 
is then applied systematically in Sees. 0] and [5] to the cases classified in Sec. El The possibility 

lr The existence of a (possibly infinite) invariant speed follows from these hypotheses and does not require an 
independent postulate. The actual value of such a speed is an experimental issue. 

2 There is no fundamental reason why it should be so [I], but there is excellent experimental evidence that 
any difference is very small. This justifies our identification. 
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for giving these theories a geometrical formulation in a four-dimensional spacetime is discussed 
in Sec. [61 Section [7] contains a few additional comments. 

We restrict ourselves to considering motion along one space dimension. The extension to 
three space dimensions is left for future investigations; see also Ref. [2]. 



2 Anisotropic relativistic kinematics 

We begin developing systematically the kinematics compatible with the principle of relativity. 
In Sec. 12.11 we limit ourselves to derive and classify the transformation laws (see also Ref. [3]), 
leaving all comments and remarks to Sec. 12.21 



2.1 Transformation law 

Consider a reference frame K, that moves with constant velocity with respect to another refer- 
ence frame /C, and assume that the corresponding temporal and spatial coordinated (t, x) and 
(i, x) are chosen in such a way that t — x — iff t = x — 0. Then, the most general coordinate 
transformation between the systems K and K, that is compatible with temporal and spatial 
homogeneity is 0, E] 

t = A(v)t + B{v)x \ . . 

x = C{v)t + D{v)x J ' K ' 

where A, B, C and D are arbitrary functions of v, the velocity!^ of /C with respect to /C. The 
parameter v belongs to some open interval J C R, containing 0. Hereafter, we shall write 
J =: ] — c_, c + [ , with c_ and c + positive constants (possibly infinite). 

The condition that x = iff x = vt, implies C(v) = A(v)v, so we can rewrite Eq. (12. ip as 

(I)= a m(!)> < 22 > 

where A(v) is the matrix 

A <"> : = ^> ( I fw) ) • < 2 - 3 > 

and we have introduced the ratios £(i>) := B{v)/A(y) and 77(f) := D{v)/A{v). Obviously, A(0) 
must be the identity matrix, so A(0) = rj(0) = 1 and ^(0) = 0. 

It has been argued by several authors that the relativity principle can be implemented at 
the kinematical level, by requiring that the transformation (12.21) be associative. Let us consider 

3 It is perhaps worth stating explicitly that the "coordinates" used in this paper (and in almost all the 
literature on special relativity) correspond to well-defined readings of time and distance, so they are not just 
arbitrary labels, but have a clear operational meaning [5]. Thus, all statements about the vague notions of 
"time" and "space" can be unambiguously interpreted in terms of the behaviour of clocks and rulers (or other 
physical systems used in measurement protocols). 

In order for the notion of velocity to make sense operationally, it is obvious that some prescription must 
have been adopted for the synchronisation of clocks in a reference frame. We do not need to specify what the 
prescription is — we only require it to be compatible with the relativity principle. 
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also a third reference frame /C, which moves with constant velocity u with respect to /C. The 
velocity u of /C with respect to K, will be given by some composition law^l 

u = $(u,v). (2.4) 

The transformation between K, and K, is 

1 X] m ( I ) . (2-5) 



so the transformation between /C and /C is 



* ) = A(«)-'A W ( \ ) . (2.6) 

By the relativity principle, this transformation must be of the same type as those that link K, 
and K to /C, so it must be possible to replace A(t>) _1 A(w) by A(u), which amounts to 

A(v)A(u) = A($(u,u)) . (2.7) 

As we shall see in a moment, this requirement allows one to determine the functions A, £ 
and rj in terms oi§ A'(0), £'(0), and r/'(0). Note that, defining as v* the velocity of a frame 
with respect to /C such that <&(v*,v) = 0, and remembering that A(0) is the identity matrix 
/, Eq. ( 12. 7ft implies A(u)A(u*) = J. Hence, imposing the relativity principle is tantamount to 
giving group structure to the set of coordinate transformations. 

Differentiating Eq. (12.71) with respect to v and evaluating the result for v = we find 

A'(0)A(w) = <f(u)A'(u) , (2.8) 

where we have defined the function 

d&(u, v 



dv 



(2.9) 



v=0 

Equation (12. 8p is a first-order, linear differential equation for A, whose solution is 

A(u) = exp (/i(V)A'(0)) , (2.10) 
where we have used the condition A(0) = /, and 

r du' 

h(u) := / ^- . (2.11) 

Jo <p{u') 



5 By the relativity principle, the velocities will all belong to the same open interval JCR, independently of 
the reference frame in which they are measured, so $ : J x J — > J. 

6 Hereafter, a prime will denote the derivative of a function with respect to its argument, with only a few 
obvious exceptions — see, e.g., Eq. (|2.11[) below. 
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The structure of A'(0) is obtained directly from Eq. (12.31) . One finds 

A'(0) = -I + M , (2.12) 

K 

where 

M:=("f ; /2 ), (2,3) 

and a := £'(()), b := r]'(0), 1/k := A'(0) + b/2. In particular, note that k is, dimensionally, a 
velocity; it can have either sign and can also be infinite. 

The matrix M has the following remarkable properties: For any natural number k, 

M 2k = (a + b 2 /A) k I , M 2k+1 = (a + b 2 /A) k M . (2.14) 

Replacing these into Eq. (I2.10p we find 

A(u) = e*")/" ( - h ( S2{u) ' 2 , , a lt\ Wo V (2.15) 
V s 2 (u) si{u) + bs 2 {u)/2 I 



where: 



+oo 1 

s iH : =E(^yj( a + &2 /4)^H 2fc ; (2-16) 



k=0 



+oo 

:= E T^TW ( a + &V4) fe M^) 2fc+1 • (2-17) 

Interestingly, s\ and s 2 satisfy the identity 

s^u) 2 - (a + b 2 /4)s 2 (u) 2 = 1 , (2.18) 

so A(tt) is, by Eq. (12.151) . just equal to e h<yU ^ K times a unimodular matrix. 
Comparing Eqs. (12.31) and (12. 15[) we find: 

A(u) = e h(u)/K - bs 2 {u)/2) ; (2.19) 

^f^-JV 1 ; (2.20) 



s 2 (u) 2, 

((m) = ow; (2.21) 

r?(u) = l + &«. (2.22) 

Note that Eq. (I2.20p contains implicitly the link between u and h(u), because its right-hand 
side depends on u only through h(u). 

Replacing Eq. ( 12.151) and the corresponding expressions for A(v) and A($(m, v)) into Eq. (12. 7ft 
we find: 

h($(u, v)) = h{u) + h{v) ; (2.23) 
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v)) = Sl (u) Sl (v) + (a + b 2 /A) s 2 {u)s 2 {v) ; (2.24) 

s 2 ($(u, u)) = si(«)s 2 (v) + s 2 (u) Sl (v) . (2.25) 

Equation (12.231) tells us that U := h(u) is additive — the canonical parametrisation for the 
group Pj. Equations (12.241) and (12.251) allow one to find the composition law between u and v, 
using repeatedly Eq. (12.201) : 

Si($(u,v)) b\ 1 w + t; + &«ti 



v ) = ' " - - = — : . (2.26) 

v ' ; \s 2 (<S>(u,v)) 2) 1 + auv K ' 

This corresponds to 

ip(u) = -au 2 + bu + 1 . (2.27) 

Finally, using the relationship s 2 (u) = si(u)/(p(u), which follows from Eqs. (12. lip . (12.161) and 
(ETTjl . together with Eqs. fl2^0l) and fl2T27D . we obtain 



s 2 («) \u 2 J —au 2 + bu + l u 2(p(u) 

This can be immediately integrated to obtain, in a neighbourhood of the origin, s 2 (u) = 
uy?(-u) _1 / 2 ; hence 

A(u)=e h ^ /K if(u)- 1/2 . (2.29) 

An alternative derivation of all these results is presented in Appendix |A] 

Equation (12.291) provides one with a criterion for finding the limit velocities — c_ and c + . 
Since f(0) = 1, it follows that J is the largest interval containing for which <p(u) > 0. The 
finite values of the limit velocities can then be found solving the equation <p(u) =0. It is not 
difficult to check explicitly that if u, v £ J, then also v) £ J, as expected for consistency. 

All these results are the most general ones compatible only with homogeneity and with 
the relativity principle. As already pointed out by Lalan [3] and others, however, not all the 
possible values for a and b are suitable for describing possible kinematics. Indeed, imposing also 
a natural pre-causality condition, that two events happening at the same place in a reference 
frame must be in the same causal relationship in any other frame, it follows that dt/dt > 0, 
i.e., that A(u) > 0, Vw £ J. This is equivalent to 

si(u) - bs 2 (u)/2 > , MueJ. (2.30) 

In order to proceed further and see what kind of constraints the inequality (12.301) imposes 
on a and b, we need to enter a tedious case-by-case analysis. Before doing so, however, it 
is convenient to change notation and use, in place of a and b, the three parameters c > 0, 
cr £ {1,-1,0}, and s £ R defined through the relations 

a = (a - e 2 ) /c 2 , 6 = 2 e/c , (2.31) 

such that a + 6 2 /4 = a jc 2 . Expressed in terms of these new parameters, the limit velocities are 

and — . (2.32) 



Of course, no limit velocities exist (not even with infinite value) for a = — 1. 
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2.1.1 Case a = 1, c < +oo 

We have 

si(v) — cosh (h(v)/c) , S2(v) = c sinh (h(v)/c) , (2.33) 

so the pre-causality condition f)2.30p is satisfied only for \e\ < 1 (which is equivalent to a > 0)0 
For |e| < 1, the equation <^(v) = has two distinct roots c + and — c_ of opposite sign, with 
c± = c / (1 =F e), and Eq. f[272Uj) gives 



c In 



l + (l + e)v/c\ 1/2 



1 - (1 -e)v/c / 

so ft, maps J = ] — c_, c + [ onto R. The transformation is therefore 



(2.34) 



c/2k 



t + (1 - £ 2 )w/c 2 



1 + (1 + e)?;/c 
1 - (1 - e) u/c,/ y/l + 2ev/c- (l-£ 2 )v 2 /c 2 
l + (l + e)w/cY /2K (1 + 2ew/c)x + t;t 



l-(l-e)v/cj y/l + 2ev/c-{l-e 2 )v 2 /c 2 



(2.35) 



obviously, a generalisation of the standard Lorentz transformation. 

In the cases e = ±1 (a = 0) there is only one finite limit velocity, equal to —1/6 = —ec/2 = 
—eC, where C := c/2 denotes the limit speed. When e = —1, J = ] — oo,C[, whereas for 
e = 1, J = ] - C, +oo [ . Equation f[2T2Q|) gives 



=eCln(l+ei;/C) 



^2.36) 



Of course, these can be considered as limit situations within the previous case, with either 
c_ — > +oo (b < 0, e = —1) or c + — > +oo (6 > 0, e = 1). 



2.1.2 Case a = 1, c = +oo (or a = £ = 0) 

The limit velocities are both infinite, so J = 1 and h(v) = v. The functions si and S2 are 
trivial: S\(v) = 1 and S2{v) = v. The transformation is 

t = e v/K t \ 

x = e v / K {x + vt) J ' (2 ' 37) 

which, apart from the factor e u / re , coincides with the Galilei transformation. These results can 
all be recovered as the limit c — > +oo of those obtained in Sec. 12.1.11 

7 Also, for e > 1 the limit velocities in (|2.32j) are both negative, and for e < — 1 they are both positive. This 
violates the condition € J. 
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2.1.3 Case o = — 1, c < +00 

As already pointed out, there are no limit "velocities". The functions s± and s 2 are 

si(v) = cos (h(v)/c) , S2(v) = c sin (h(v)/c) , (2.38) 

and 



v/c 



tan (h(v)/c) 



(2.39) 



1 + ev/c 

The transformation is better expressed in terms of the canonical parameter V = h(v), since the 
use of v would require some sign ambiguities in order to describe it completely: 



•t = e v/K 



x = e 



V/k 



cos (V/c) - e sin (V/c) j ct - (1 + e 2 ) sin (V/c) x 
sin (V/c) ct + ( cos (V/c) + £ sin (V/c) ) x 



(2.40) 



This generalises rotations in a Euclidean plane. Pre-causality is violated for all choices of e. 



2.1.4 Case o = 0, c < +00, e 7^ 
The function ip is 

cp(v) = (l + bv/2) 2 , 
and there is only one limit "velocity", C := —2/6. One finds 

s x (t)) = 1 , s 2 (v) = h(v) 



1 -v/C 



(2.41) 



(2.42) 



and the transformation is 

Qh(y)/n 

t 



1 -v/C 



3 ft(«)/ 



K (1-2 v/C) x + wt 



x = e wc^ = e 



y/K ((l- V/c) x + Vt) 



^2.43) 



Note that all these expressions can also be obtained by taking the limits c — > +00, £ — > ±cxd, 
while keeping the ratio e/c finite, of those found in Sec. 12.1.31 Pre-causality is violated. 



2.2 Comments 

Summarising, the only cases that correspond to physically acceptable kinematics are those 
considered in Sees. 12.1.11 and 12.1.21 We now comment on the properties of the generalised 
composition law for velocities (Sec. 12.2.11) . on the physical interpretation of the parameters e 
and k (Sec. 12.2.21) . and on the extension to three spatial dimensions (Sec. I2.2.3j) . 
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2.2.1 Velocity composition law 

The function <3> given in Eq. (12.261) satisfies the following properties: 

$(u,0) = $(0,u) = u, VmG J; (2.44) 

VuG J, 3u*eJ such that $(«*, u) = u*) = ; (2.45) 

$($(«, i;),u>) = $(it, , Vu,v,w G J ; (2.46) 

$(«, v ) = $(v , «) , Wu,v E J . (2.47) 

Hence, on writing u © w := $(w, t>), \/u,v G J, Eq. (I2.26P defines the composition law of an 
Abelian group (J, ©), with neutral element and inverse u* of a generic element u G J defined 
by Eq. (!2.45p Fl This can also be derived straightforwardly from simple kinematical arguments 
and the principle of relativity [lj. Remarkably, however, the associative and commutative 
properties (I2.46P and (I2.47P do not hold in general for the composition law of velocities along 
arbitrary directions in more than one spatial dimension [7]. 

Since the function $ is infinitely different iable in both its arguments, (J, ©) is a C°° one- 
dimensional Lie group. Any smooth, connected one- dimensional manifold is diffeomorphic 
either to R, or to one of the real intervals [0, 1] and [0, 1[ , or to the circle S 1 (see, e.g., Ref. [8]). 
Hence, every one-dimensional Lie group is isomorphic either to (R, +) or to [0, 1[ with the 
addition modulo 1, according to its connectivity. This implies that, in our case, there exists 
a C7°° additive function defined on J, taking values either in R or in [0, 1[, for which is a 
fixed point0 Of course, this coincides with the function h previously defined. Since h is one-to- 
one, h(v) (called rapidity in the literature on special relativity) can be taken as an alternative 
mathematical representation for the physical notion of velocity, instead of the more common v. 
This idea is supported by the observation that h(v) also admits an operational definition [10J. 

Suppose that a point has velocity u with respect to the reference frame /C. Its velocity u 
with respect to K, is also given by the composition law ( 12.261) . as one derives straightforwardly 
from the transformation equations (12.11) . with and n(v) given by Eqs. (12.211) and (12.221) . 
respectively, with the obvious replacement u — > t>0 However, now only the velocity v (of 
/C with respect to /C) is forced to belong to the interval J, because the moving point is not 
necessarily associated with some material particlef"1 

The function (p contains all the information needed to specify $. Its meaning can be found 
by expanding u to the first order in v. 

u = u + (f{u)v + 0(v 2 ) . (2.48) 

This is the composition law between an arbitrary velocity u and a velocity v with small mag- 
nitude. Since Eq. (12.441) implies y?(0) = 1, at very small speeds one always recovers Galilean 



8 In general, u* — —u/ (1 + bu), so u* — —u only when 6 = 0. 

9 See Ref. [9] for an alternative proof of this theorem, with applications to relativistic kinematics. 
10 Note that the function A does not enter the velocity composition law. 

1 It could correspond to a purely geometrical occurrence like, for example, the intersection between two 
moving straight lines, which can have an arbitrarily high speed even if the lines move rather slowly, provided 
the angle they form is small enough. 
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kinematics. However, unless <p(u) = 1, deviations from the Galilean composition law can al- 
ways be detected if one measures the speed of a fast-moving object, even for small values of 
the relative velocity v between frames. 

Another interesting property of if is that if some finite velocity, say C, is invariant, then 
<f(C) = 0. This follows immediately by applying Eq. ( 12. 9ft to the condition 

$(C» = C, VveJ, (2.49) 

which expresses the invariance of C. Hence, the limit velocities for reference frames coincide 
with the invariant velocities^ 

2.2.2 Conventional and real anisotropy 

For no n- vanishing values of the quantity e e] — 1, 1[, the invariant speeds c± along the two 
orientations of one-dimensional space differ from each other [T7] . On the other hand, the average 
(two-ways) speed of a signal travelling at the invariant velocity along a round trip turns out to 
be equal to the e-independent parameter c. For e — > ±1, one of the two invariant speeds tends 
to infinity, while the Galilean composition law is recovered in the limit c — > +oo. A further 
requirement of spatial isotropy (or better, its one-dimensional counterpart - - the physical 
equivalence of the two orientations in the one-dimensional space) enforces, not surprisingly, 
the equality between c+ and c_ [3]. Hence, anisotropy seems to be related to the possibility 
of having c + ^ c_. However, whether the one-way speed of light is a physically meaningful 
quantity or merely a conventional one, is the matter of a long-standing debate [TTJ [T8] . 

This issue is inextricably linked to another one, concerning the conventionality of clock 
synchronisation [T3"j [TO] [T8] . Any measurement of the one-way speed of light requires, in order 
to be performed, a prior synchronisation of distant clocks. Vice versa, one may argue that any 
synchronisation procedure is equivalent to a stipulation about the value of the one-way speed 
of light @ This creates a circularity which does not allow any escape. Indeed, as it is evident 
from the first of equations (12.351) . e also encodes the convention adopted for synchronising 
clocks. The value e = corresponds to the usual (Einstein's) procedure, while other values 
account for Reichenbach's generalised synchronisations [T3J [TBI OS] • Thus, for synchronisation 
procedures different from Einstein's, c + 7^ c_, whereas the value of c does not depend on the 
synchronisation procedure adopted. 

At a first sight, this circumstance appears puzzling. Isotropy, being related with the equality 
of the one-way invariant speeds c±, can be regarded as a consequence of Einstein's synchroni- 
sation procedure, and it has been argued convincingly that the choice of such a procedure is 
merely conventional, in spite of some claims to the contrary^ On the other hand, it seems 
reasonable to believe that isotropy (or anisotropy) is a physical property, which cannot be 
implemented or altered just by a stipulation. 

12 Note that the possibility for the existence of invariant speeds has been derived as a kinematical possibility 
only from the postulates of relativity, homogeneity, and pre-causality. This approach to relativistic kinematics 
was pioneered by von Ignatowsky in 1910 [llj . and was later rediscovered many times in different ways 3,9, 
[T2][T3]. See also [14] for a rigorous treatment, and [HI [16] for clear presentations at a textbook level. 

13 These statements are less tautological than it may seem at first, as one can synchronise without using light. 

14 See PS US EH] for the debate. 
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The resolution of this conundrum is that e = is only a necessary, but not sufficient, 
condition for isotropy. The transformation (12.351) contains also the parameter k, and if |ac| < 
+oo one can experimentally distinguish between the two spatial orientations even if e = 0. 
For example, if one considers two clocks moving at the same speed along opposite directions, 
these clocks will not delay by the same amount, if \k\ < +oo. Hence, anisotropy can manifest 
itself through physical effects, independent of the choice of synchronisation. One should thus 
distinguish between a real, physical anisotropy, measured by the parameter k, and a purely 
conventional one — that can be introduced or gauged away simply by a stipulation — measured 
by e. From this point of view, the fact that one may not use Einstein's synchronisation is a 
trivial one, and has no physical content. On the contrary, the fact that one can adopt it in any 
inertial frame is physically non-trivial. 

An analogy to this situation is provided by the principle of inertia. The statement that, in 
an inertial frame, a force-free particle moves along a straight line at a constant speed contains 
two elements, very different in nature. One, that the motion takes place along a straight line, 
is a physically testable prediction, since the notion of a straight line is well defined in the 
Euclidean geometry that one presupposes valid when formulating the principle. The other, 
that such motion is uniform, is a matter of convention. Of course, one could choose the "time" 
variable in such a way that the motion is not uniform (the analog of choosing a synchronisation 
different from Einstein's), but this generalisation will not lead to new physical phenomena - 
only to a horrendous complication in the formulation of the laws of mechanics. Again, the 
relevant fact is not that one can make an absurdly complicated choice of time, but rather 
that one can make a choice that simplifies life. Since a change in the time variable does not 
entail new phenomena, we are confronted with a mere gauge, and the wisest choice is to use 
the simplest possible gauge. As Misner, Thorne and Wheeler concisely and effectively wrote: 
"Time is defined so that motion looks simple" [20] . We could paraphrase them saying: "Clocks 
are synchronised so that physics looks simple." 

Intriguingly, there is a possibility that the value of k also reflects a mere convention. The 
generalised Lorentz transformation (12.351) with e = is a standard Lorentz transformation 
accompanied by a global dilatation by a ^-dependent factor, dilatation which can be attributed 
to a rescaling of the time and distance units in the reference frame /C. Such a rescaling could be 
an effect of anisotropy on clocks and rulers, but could also be induced artificially, by a suitable 
choice of units in different frames. In any case, it can be compensated through a rescaling of 
units in the reference frame /C0 

There is a way out of this difficulty, though. The anisotropic scale factor in Eq. (I2.35j) 
is very special, and requires a carefully tailored choice of units in fC in order to be produced 
artificially. In common experimental practice, there is a well-defined (although not explicitly 
stated) procedure for building units in a frame /C, say. They are either constructed directly 
in K, following some standard instructions; or they are boosted from another frame /C, where 
they have been built using the standard instructions. (That these two operations produce the 
same units in /C, is the so-called principle of the boostability of units [IS].) Thus, although the 

15 This operation does not destroy the group structure, and is therefore compatible with the principle of 
relativity. On the contrary, eliminating the Lorentz factor by a rescaling of units would not preserve the group 
structure, hence would imply a violation of the principle. 
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choice of units in different frames is, in principle, free, in practice it is always made in the same 
way — a way, moreover, which does not seem capable of introducing an anisotropy. Hence, 
any detection of the anisotropy allowed by Eq. (I2.35P will be regarded as a real physical effect, 
rather than as an artifact of conventions! 16 ! Such an effect could be gauged away only at the 
price of changing everyone's measurement habits. 



2.2.3 Three-dimensional case 

Although our treatment was restricted to the case of one spatial dimension, it is easy to extend 
it to a three-dimensional space. We do not provide a systematic generalisation here, but 
nevertheless we wish to present a physical argument that shows what one could expect. (See [2T1 
[22] for more details.) For simplicity, we choose Einstein's synchronisation (e = 0). Consider 
a light-beam clock, of the type that is commonly used in pedagogical treatments of special 
relativity (see, e.g., Ref. [23]). Basically, such a clock is made of two mirrors, one along the 
x-axis, the other displaced from it along the perpendicular direction. Let us denote by y and y 
the distance between the mirrors in the reference frames K, and /C, respectively. The time taken 
by light to make a complete two-way trip between the mirrors is, in /C, simply At = 2y/c. The 
corresponding time according to /C is then given by the first equation in (I2.35P with 5 = 0: 

\1-V/Cj y/l - V 2 /C 2 

But At can also be found directly from the equation 



At = - \ (vAt/2) 2 + y 2 , (2.51) 
c v 

which gives 

Ai= , 2V/ ° =. (2.52) 

y/l ~ V 2 /C 2 

Combining Eq. fl2"3UD with At = 2y/c and Eq. f[2"32|) . we find 



1 — v/c / 

instead of the usual y = y. Thus, anisotropy introduces also a transformation for distances 
along transverse directions. 



3 Foundations for dynamics 

We now present a method for finding the basic dynamical quantities (Lagrangian, energy and 
momentum) for a free particle, that are compatible with the kinematics presented in Sec. [2j 

16 The situation is different for the transformation (|2.40p . When s = this is a rotation, anticlockwise by an 
angle 8 := V/c, of orthogonal axes in a Euclidean plane, accompanied by a global dilatation by the factor q 6c / k . 
A true anisotropy cannot produce effects for 9 = 2it, so the factor e ec / K can only be due to a devious choice of 
units in the frame K-. 
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3.1 Lagrangian 

The Lagrangian L for a free particle in an inertial frame can be determined, in Newtonian 
dynamics, by the following argument |24J. Space and time homogeneity requires that L do not 
depend on the particle position and on time, so it must be only a function of velocity. By the 
relativity principle, we have that in the inertial frames K, and K, this function must be the same, 
so we shall write L(u) and L(u), respectively, where u — u+v. However, these two Lagrangians 
must lead to the same equation of motion, so they can differ by a total derivative with respect 
to time of a function / of the particle coordinate and of time. Hence 

T ( . x t / \ I df(x,t;v) df(x,t;v) . . 

L(u + v)=L(u) + — u + . (3.1) 

Since the left-hand side does not depend on x and t, the same must happen for the right-hand 
side, so the quantities df/dx and df/dt can actually be functions of v only. Denoting them as 
a(v) and /3(v), we have therefore 

L(u + v) = L{u) + a(v) u + (3(v) . (3.2) 

Taking the derivative of Eq. (13.21) with respect to v, and evaluating the result for v = 0, we get 

L'(u) =a'(0)u + (3'(0) . (3.3) 

Integrating, one finds that L(u) is a polynomial of second degree in u. Discarding the part 
linear in u, which does not contribute to the equation of motion, one ends up with the usual 
result: L(u) is proportional to u 2 . Note that this implies isotropy of space, although no such 
hypothesis has been explicitly used. 

It is very easy to see that one cannot run this argument in a straightforward manner to 
get the Lagrangian corresponding to a generic composition law (12.41) . Indeed, even for the 
composition law of special relativity one would not recover the standard expression for the 
Lagrangian. The reason lies in the fact that Eq. (13.11) is appropriate only for a theory in which 
the time coordinates t and t in K and /C, respectively, coincide. As we saw in Sec. El this 
happens only in isotropic Galilean kinematics. 

This difficulty can be overcome by working in the extended configuration space of the 
particle, where position and time are both treated as Lagrangian coordinates, evolving in terms 
of a parameter 9 [23]. The condition that the Lagrangian, expressed in the inertial frames JC 
and /C, lead to the same equation of motion is then 

dt , / dx , di\ dt r f dx , dt\ , , dx „. . dt , 4 . 

U ' W = dfl* U 1 dfl) + a(v) Te + «"> To • < 3 ' 4 > 

where x and x denote the particle position in K, and /C, respectively. Equation (13. 4p can be 
rewritten as 

g(u, v)L($(u, v)) = L(u) + a(v) u + (3(v) , (3.5) 
where we have defined the function 
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such that g(u, 0) = 1. 

We can now proceed as in the Galilean case, by taking the derivative of Eq. (13.51) and 
evaluating it for v = 0. The result is the following first-order ordinary linear differential 
equation for the function L(u): 

<p(u) L'{u) + p{u) L(u) = a'(0) u + /3'(0) , (3.7) 

where 



dg(u,v) 

PW ■ = 



9V v=Q 

Using Eqs. Q, Q, flCTj) and (jMH) into Eq. (j3TB]) . we now find immediately 



g(u,v) = (l + auv)A(v) . (3.9) 

Hence, remembering that A(0) = 1, 

p(u) = au + A'(0) = au + 1/k- b/2 , (3.10) 

and the differential equation ( 13. 7ft becomes then 

(l + 6u-aM 2 )L'(M) + (aM + l/K-6/2)L(M) = a'(0) u + /3'(0) . (3.11) 

It would now be straightforward to find the Lagrangian corresponding to the various possible 
kinematics classified in Sec. |2j However, we find it interesting to show that the same functions 
ip and p can be obtained without any explicit reference to the derivation in Sec. 12.11 simply 
requiring that the dynamics based on the Lagrangian L be compatible with the existence of 
elastic collisions between asymptotically free particles. 



3.2 Momentum and energy 

Given the Lagrangian L(u), the particle momentum and energy^! are easily found by the 
standard relations: 

p(u) = L\u) ■ (3.12) 

E{u) = up(u) - L{u) . (3.13) 

Noteworthy, expressions for p(u) and E(u) can also be derived following the procedure 
in Ref. which relies directly on the existence of elastic collisions and on the relativity 
principle Jlj instead of starting from a Lagrangian. Let E(u) be the energy of a particle with 
velocity u in an inertial frame /C. In an inertial frame, E(u) is conserved for a free particle, 
because u is constant, by the principle of inertia. Then, the total energy — defined as the sum 
of the energies for the individual particles — is conserved also for a system of noninteracting 
particles. 

17 Hereafter, whenever we refer to "energy" we mean the sum of kinetic energy and a possible rest energy. 
18 Basically, a straightforward generalisation of an argument originally due to Huygens [26]. See also [27l [28] 
for similar developments. 
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Let us assume that there are spatially localised interactions between particles which do not 
change the total energy. Then, during such an interaction between two particles with initial 
velocities u[ , u 2 9 , and final velocities uf , u 2 ' 

E^uf) + E 2 {uf) = E x {uf) + E 2 (uf) . (3.14) 

(Of course, the energy may depend on some invariant parameter characterising the particle, in 
addition to its velocity. For example, in Newtonian dynamics it depends on the particle mass. 
We keep track of this dependence through the indices 1 and 2 on £.) With respect to another 
inertial frame /C, 

E^uf) + E 2 {uf) = BktfP) + E 2 (4 ] ) , (3-15) 

where the various and are given by the composition law for velocities (12.41) . and we have 
used the same functions E\ and E 2 in both reference frames because of the relativity principle. 
We now expand the generic functions in Eq. (13 . 1 5[) around v = 0, and use Eq. (12.481) to get 

E(u) = E{u) + E\u) <p{u) v + 0(v 2 ) . (3.16) 

Doing this for each term in Eq. (13.151) and using Eq. (13.141) . then dividing by v and taking the 
limit for v — > 0, we obtain 

E[(uf) ip{uf) + E' 2 {uf) ip{uf) = E[(u { P) v(u?) + E' 2 {uf) tp{uf) . (3.17) 

Hence, there is another additive quantity which is conserved, in addition to energy. For a single 
particle, the most general expression for such a quantity is 

p(u) = A tp(u) E'{u) + /i E(u) + v , (3.18) 

where A, fi, and v are quantities independent of u. The function p(u) given by Eq. (I3.18P 
with A = l,/i = z/ = coincides with linear momentum!^! both in Newtonian and Einstein 
mechanics PQ, and we shall retain such an interpretation for the more general expression above F°l 

3.3 Compatibility 

If one assumes that p and E can be derived from a Lagrangian according to Eqs. (13.121) and 
(13.131) . Eq. (I3.18P can be converted into a second-order differential equation for L: 

Xu(f{u)L"{u) + {fiu-l) L'(u) - /iL(u)+v = . (3.19) 

In order for the treatments in Sees. I3.l l and !3.2l to be mutually compatible, Eqs. (13.71) and (I3.19P 
must have the same content. To compare them, let us first rewrite Eq. (13. 7p as a second-order 
differential equation. By taking a derivative with respect to u, then multiplying by u, we obtain 

u (p(u) L"(u) + u (ip'(u) + p(u)) L'(u) + u p'{u) L{u) - a'(0) u = . (3.20) 

19 Note that with this identification, linear momentum turns out to be (correctly) a one-form rather than a 
vector [2]. 

20 No new conservation laws can arise at higher orders in v, as argued by Levy-Leblond [28] . See also Ref. pQ for 
an explicit proof when E is only a function of u 2 (which, however, is not the case in the presence of anisotropy) . 
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We can now replace the last term in Eq. (I3.20p . o/(0)u, by using again Eq. (13. 7p . The final 
equation is 

u (p(u) L"(u) + {u cp'{u) - <p(u) + u p{u)) L'(u) + {u p'{u) - p{u)) L(u) + 0'(O) = . (3.21) 

Equations H3 . 191) and (13.211) coincide if the following relations hold: 

/3'(0) = u/X ; (3.22) 

up'(u) - p{u) = -At/A ; (3.23) 
u (f'(u) — (f(u) + u p{u) = pu/\ — 1/A . (3.24) 
The most general form of the functions p and ip allowed by Eqs. ( 13.231) and ( 13.241) are 

p{u)=au + p/\ (3.25) 

and 

tp(u) = -a u 2 + b u + 1/A , (3.26) 

with a and b arbitrary constants!^] Moreover, Eq. (I2.44p implies ^(0) = 1, so A = 1. Hence, 
the treatments in Sees. 13.11 and 13.21 are compatible only if: 

p(u) = au + p; (3.27) 

v (u) = - au 2 + bu + l- (3.28) 

v = (3'(0) . (3.29) 

Remarkably, the structure of p and ip emerges by the requirement that the treatments in 
Sees. 13.11 and 13.21 be compatible, with no independent considerations about kinematics. Note 
that the expression for <p(u) coincides with the one given by Eq. (12.271) . obtained on purely 
kinematical grounds. Comparing now Eq. (13.271) with Eq. (13. lOj) . we are led to the identification 

p = A'(0) = 1/k - b/2 . (3.30) 
3.4 Mass, rest energy, rest momentum 

Evaluating Eq. (13. lip and its first derivative at u — 0, we can express the quantities ct'(0) 
and /3'(0) in terms of a, b, p, and of the three parameters m := L"(0), p := p(0) = £'(0), 
and E := E(0) = — L(0), that represent the particle's mass and possible rest momentum and 
energy, respectively!^! The results are 

a'(0) = m+ (b + p)p - aE (3.31) 



21 The constants a and b have, a priori, nothing to do with those introduced in Sec. 12. 1[ but we shall soon 
discover that they actually coincide. This justifies using the same letters in the notation. 

22 One could distinguish between the momentum p and a "kinetic momentum" p — pq, just as one usually 
distinguishes between the energy E and the "kinetic energy" E — Eq. 
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and 

/3 , (p)=p -fiE , (3.32) 
so the differential equation for L is, finally: 

(l + b u — a u 2 ) L'{u) + (au + fi) L{u) = (m + {b + u) p — a E ) u + p — fi E . (3.33) 

This is the equation we shall use as a starting point in the next section. 

It is important to note that the parameter m coincides with the Newtonian mass, thus 
justifying our setting m := L"(0)o This can be seen by expanding the Lagrangian and keeping 
only the leading order terms in m/c and u/k. The result is 

L(u) = ^mu 2 +p u- E + O(u 3 ) , (3.34) 

which to the first significant order in u describes indeed a Newtonian particle with mass m. 

Instead of solving Eq. (13.331) . one could proceed as follows. Differentiating Eq. (13.131) . and 
using Eq. (I3.12p . one obtains the well-known relation 

dE(u) =udp(u) , (3.35) 

expressing the fact that the change in the energy of a particle equals the work done on it 
- an elementary property that holds not only in Newtonian and Einstein mechanics, but in 
any dynamics that admits a Lagrangian formulation. Combining Eq. (13.181) with A = 1 and 
Eq. (13.351) . one obtains a single differential equation, which can be solved to find the expression 
for E(u). Inserting the latter into Eq. (13.181) . one finds also the expression for p(u). 

Finally, let us notice that momentum and energy must be proportional to the particle mass 
m. This follows from the requirement that momentum, energy, and mass are all additive 
quantities. Denoting by p(m, u) and E(m, u) the momentum and energy of a particle with 
mass m and velocity u, we have thus, for a system of two particles with masses mi and m<i that 
move with the same velocity u, 

p(mi,u) + p{m 2 , u) = p(m x + m 2 , u) , (3.36) 

and 

E(mi, u) + E(rri2, u) = E{mi + m 2 , u) , (3.37) 

from which the said proportionality easily follows. As a byproduct, the rest momentum p and 
the rest energy E Q turn out to be also proportional to m. 

4 Anisotropic relativistic dynamics 

We now write the explicit expressions for the basic quantities that appear in an anisotropic 
relativistic dynamics. Although the considerations in Sec. 12.2.21 strongly suggest to set e to 



23 It might have been logically possible that L"(0) were equal to the Newtonian mass multiplied by a function 
of c and k that reduces to 1 when both these parameters tend to infinity. 
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zero, and experimental evidence suggests k > c, in the following we shall consider all the values 
of the parameters that correspond to possible kinematics, in order to remain as general as 
possible. The differential equation (I3.33j) has two different solutions when \k\ ^ c and when 
\k\ = c, so we must consider these two cases separately. 

4.1 Case \k\ ^ c 

The solution of Eq. (13 .33|) that satisfies the condition L"(0) = m is 

L(u) = - 1/9 m n/ 9 e- fc <">/« <p(u)^ + ( Po -m }' R - e '\ \ u-( Eo 



1/c 2 -1/k 2 rv ' V 1/c 2 -1/k 2 J \ J 1/c 2 -1/k 2 

'(4.1) 

The expressions for momentum and energy can then be obtained by Eqs. (13.121) and (13.131) : 

*' = l 7 ^ e " M " ,/ >("»" 1/2 (^" + ^j)K^- m T7^) ; (42) 



(4.3) 



4.1.1 Anisotropic Einstein's dynamics 

Replacing Eq. <^Mj into Eqs. (jPjl- fO) we find: 

L(u)= - ~rTT~~rT^> ( !:'!:!? (1 - (1 - ») »/cr (1 + (1 + ») 



1/c 2 - 1/k 2 VI + (l+e)«/c 
/ 1/k — e/c \ /_ m , 

+ ("°-"' 1 / tf - i ; K 2 )»-K- 1/c2 _ 1/K J ; en 



m 

p(u) = 



1 - (1 - e) u/c\ c/2k (1 - e 2 ) m/c 2 - e/c + 1/k 



1/c 2 - 1/k 2 V 1 + (1 + e) u/cj (1 - (1 - e) u/cf 2 (1 + (1 + e) u/c) 



1/2 



1/k - e/c , t . 

+ P0 ~ m 1/C 2 - 1/K 2 ; (45) 



E u) = 



l-{l-e)u/c\ c/2K 1 + (1/k + e/c) u 



1/c 2 - 1/k 2 Vl + (1 + e) «/c/ (1 - (1 - e) m/c) 1/2 (1 + (1 + e) u/c) 



1/2 



771 

These equations generalise the basic expressions of Einstein's dynamics (corresponding to e = 0, 
\k\ = +oo) to the anisotropic case. The cases in which anisotropy is merely due to a convention 
(e 7^ 0) and in which it is an intrinsic physical property (\k\ < +oo) are both covered. Figured] 
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Figure 1: Comparison between the expressions given by Eqs. (|4.5p and (|4.6p for k = 2 c (solid lines) 
and |k| = +oo (dashed lines). In both cases, the conventional anisotropy parameter e has been set 
equal to zero, and the values of the quantities at rest have been chosen according to Eqs. (14.7)) and 
(|4.8p . The plot on the left represents momentum, the one on the right represents energy. 



shows a comparison between these expressions and those valid assuming isotropy. Note that, 
even in the presence of anisotropy, the minimum value of E is still attained for u = 0. 

Remarkably, Eqs. ( 14. 5ft and (14. 6 ft suggest both the possibility for a non-vanishing rest mo- 
mentum 

1/k — e/c 
1/c 2 - 1/k 2 

and a modification of the celebrated equation Eq = mc 2 into 



j. / n, — c / (_ . . 

P° = m TT3. 77^' ( 47 ) 



m c 2 



E ° = ■ (48 > 

Note that by Eqs. (14.71) and (14.81) the right-hand side of the differential equation (13.331) vanishes, 
so these choices for p and E are equivalent to setting the constants a'(0) and P'(0) in Eq. (13. lip 
equal to zero. Of course, Eq. ( 14.81) is physically sound only if one can argue that, for a system 
of particles, the total energy in the centre-of-momentum frame stands in the same relation with 
the mass of the system, regarded as a single unit. This is shown in Appendix iBl 

Finally, let us consider the extreme cases of strong conventional anisotropy, e = ±1. Now 
one of the two invariant speeds c+ and c_ is infinite, and it is convenient to write, as we 
already did in Sec. 12.1.11 C = c + or C = c_, according to whether b is positive or negative, 
respectively. This amounts to setting b = e/C, or equivalently c = 2 C. The expressions for 
the basic dynamical quantities can either be obtained directly from Eqs. (I4.1l) - (l4.3p . or from 
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Eqs. fl4.4|) - fl4.6p with the replacements e = ±1 and c = 2C: 

T , m ( eu\ l / 2 ~ £C / K 

L(U)= - l/AC 2 - lp V + -C) 

(Tfl \ ( 771 \ 

Po + e/2C + l/. ) U -{ E °- l/4Ci-l/^ ) ; (49) 

^ ) = -7^cTy-A 1 + ^) +p ° + 7l2cTTT K l (410) 



1/AC 2 -1/k 2 \ CJ \ \2C kJ J u 1/AC 2 -1/k 2 

4.1.2 Anisotropic Newtonian dynamics 

We now consider the case a = e = (or, alternatively, c = +oo). On replacing a = b = 0, 
fM — 1/k, h{u) = u into Eqs. fl4.1l) - fl4.3p . one finds: 

L(u) = -m k 2 (1 - e~ u/K ) + (po + m k) u - E ; (4.12) 

p(u) = m/c (1 -e~" /K ) +p ; (4.13) 
E(M) = m/t 2 (1 -e" n/K ) -mKae^ + Eo . (4.14) 

Note that for \u\ |k| one recovers the usual expressions of isotropic Newtonian dynamics 
(apart from the constants po and Eq). Figure [2] shows a comparison between these expressions 
and those valid assuming isotropy. Both p and E tend to asymptotic values as u tends to +oo 
if k > 0, or — oo if k < 0. This behaviour, although mathematically interesting, is however 
irrelevant from a physical point of view, because it concerns the high-speed regime of a theory 
which is, in fact, only a low-speed approximation. 

Equations fl4.12l) - fl4.14l) can be obtained, formally, as the limit for c — ► +oo of the corre- 
sponding expressions in anisotropic Einstein's dynamics, Eqs. fl4.4l) - fl4.6p . However, if the limit 
is taken after the choices (14.71) and (14. 8[) for Eq and po are made, one gets a rather bizarre 
behaviour for the rest energy and momentum in the Newtonian limit, as these quantities would 
diverge as |«| — > +oo. Of course, in drawing this unphysical consequence one ignores the fact 
that in the real world c is finite, and that experiments give \k\ > c, so one cannot really study 
the case c — > +oo keeping |«| finite. 

4.2 Case \k\ = c 

The exceptional cases with k = ±c have of course no physical interest; nevertheless, we present 
the results for the sake of completeness. The solution of Eq. (13.331) can be written in a unified 

way as 

r , s mc 2 ( (\ e\ \ , fl — (1/ 'k — e/c)u\ , . > , J „. 

= — ( 1 " (n -c) U ) ln { l + \l/ K J/c)u + m K/2) U - Ea • (415) 
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Figure 2: Comparison between the expressions given by Eqs. (|4.13p and (|4.14|) for < k < +oo (solid 
lines) and \ n\ = +oo (dashed lines). In both cases, the values of the quantities at rest have been chosen 
equal to zero. The plot on the left represents momentum, the one on the right represents energy. Note 
the horizontal asymptotes of p and E as u — > +oo (displayed as thin solid straight lines). 



where again Eqs. (I2.3ip and f 1 3 . 3 j) have been used, and m = L"(0). The momentum and energy 
are then, respectively: 



p(u) 



ttlk (1/k, + e/c)u mc 2 



E(u) 



2 l + (l/« + e/c)u 4 
1 



--^lln 

K C 



(l//t — e/c)u 



niKU 



mc 



In 



1 + (1/k + s/c)u 
1 — (1/k — e/c)u 



+ E . 



(4.16) 



(4.17) 



2 1 + (1/k + e/c)u 4 V 1 + (l/« + e/c)u. 

Note that these expressions cannot be obtained by taking the limits n — > ±c of the correspond- 
ing expressions in Sec. I4.1[ because such limits do not exist. 

In the extreme cases of strong conventional anisotropy (e = ±1, c = 2 C), these expressions 
reduce to: 



L(u) 



p(u) 



emnC 



- ~ ^) A In (l + ^t) + (po + m k/2) u - E ; 



1/k + e/2C)u emnC fl e \ f eu^ 
2~l + (l/n + s/2C)u + ^2~{K~2c) ^ + 



mnu 



1 



emnC 



eu 



(4.18) 
(4.19) 
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5 Hamiltonian and dispersion relation 

Comparing Eq. (13.351) with Hamilton's equation 

« = ^, (5.1) 
dp 

allows us to identify the Hamiltonian H(p) as E(u(p)), up to a velocity-independent additive 
term that we set equal to zero. (Of course, the same expression for H can be obtained as the 
Legendre transform of L.) 

In principle, one could find H directly, solving a differential equation as we did for L in 
Sec. HI Indeed, by a further differentiation of Hamilton's equation (15. ip we get 

du = ^dp. (5.2) 
dp z 

Combining this with Eq. (I3.18P for the relevant case A = 1, we arrive at the following differential 
equation for H(p): 

d 2 H dH fdH\ 

In general, however, Eq. (15.31) is not easy to solve. Therefore, it is more convenient to rely on 
the expressions for E{u) and p(u), trying to eliminate u to obtain an implicit relation between 
energy and momentum, from which the Hamiltonian can, in principle, be extracted by local 
inversion. This implicit relation is the particle version of a dispersion relation. 

5.1 Case \k\ ^ c 

5.1.1 Anisotropic Einstein's dynamics 

It is convenient to define the new quantities 

/ s / c 

p(u) := p(u) -p + m _ 2 (5.4) 

and 

T?7 

E{u) := E{u) -E 0+ — _ — , (5.5) 
1/k — e/c , 

pi - 0) = m W^W (5 ' 6) 

— JJl 

m = W^W ' (57) 

independently of the choices for po and Eq. (Hence, p and E coincide with momentum and 
energy once the natural choices for p and E have been made, as in Eqs. (14. 7p and (14.81) .) 
Then, let us form their linear combinations 

(l-e)E-pc, (l + e)E+pc } (5.8) 



which are such that 



and 
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and multiply and divide these by each other. In this way, we obtain two expressions that both 
contain u in the combination 

1 - {l-e)u/c 
l + (l+e) u/c ' 

so at the end we can form a w-independent expression: 



1— c/k / \ 1+c/k 

\-e)E -pc\ ((l + e)E + pc 



m 2 c 2 



An alternative form of this equation, which allows also for negative values of E, is 

1+c/k 



(5.9) 



1 -e)E-pc) (1 + e)E +pc 



1-c/k 



Tfl 1-c/k qI-c/k 



1 + c/ 



(1/c- 1/k) (1/c+ l/«) 



(5.10) 



This is the sought-for dispersion relation, that generalises the well-known equation E 2 —p 2 c 2 = 
m 2 c 4 to the case in which space is anisotropic. Figure [3] shows a comparison between these two 
expressions. 




Figure 3: Comparison between the set of pairs (p,E) that satisfy Eq. (|5.10j) for k = 2 c (thick solid 
line) and \k\ = +oo (dashed line). The conventional anisotropy parameter e has been set equal to 
zero, and the asymptotes are also displayed (thin solid straight lines). 



To the first order in c/k, Eq. (I5.9p leads to the approximate relation 

CI - :)E-pr) (CI ■-;) E ■->•) ^ m 2 c 4 ( 1 + - In U l ~ £ ) E -P c 

K \(l + e)E + pc 



(5.11) 
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which can be used in order to place experimental bounds on c/k. 

In the cases of strong conventional anisotropy (e = ±1), using the same notations as in the 
final part of Sec. 14.1.11 and eliminating u, we can directly write the Hamiltonian: 

^-H 1 + ^fa^)'* B )- (5 ' 12) 

Of course, this can also be recovered by Eq. (15.101) with the appropriate substitutions. 

5.1.2 Anisotropic Newtonian dynamics 

The Hamiltonian is 

H(p) = E Q + K(p-p )+mK 2 (l-^^-)ln(l-^^) 

\ mK J \ mn J 

= E + ^-(p-p ) 2 + O((p-p ) 3 ) . (5.13) 
2m v 

Figure H] shows a comparison between the cases with and without anisotropy. 



H-E 






P-Po 



Figure 4: Comparison between the set of pairs (p—po,H — Eo) that satisfy Eq. (j5.13j) for < k < +oo 
(solid line) and |k| = +oo (dashed line). Note that the curve for a finite k possesses an end-point (P 
in the diagram), which corresponds to the finite asymptotic values of p and E as u — > +oo. 
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5.2 Case \k\ = c 

By the same procedure as in the previous section we find 

1 s\ ~ mu 

} E —p = — -— 5.14 

k cj y l + (l/« + e/c)u v ' 

and 

1 e\ ~ mc 2 fl - (1/k — e/c)u\ /r . 

- + -\E + p= In ^ ——- — ] , (5.15) 

k c J y 2k \1 + (1/k + e/c)u) V ; 

where now E = E — E and p = p — Pq. Extracting u from Eq. (I5.15P and replacing it into 
Eq. (15.141) one finds, after trivial algebra, the dispersion relation: 

e\ ~ vtik vtik ( 2 k ( ( 1 e\ ~ \ \ 

__j B _ J)+ ___ exp ^_^_ + _j B+p -jj = 0. (5.16) 

For strong conventional anisotropy we find 

^-^j 5 -^ — exp (-^(^ + ^J £+!, Jj =0 ' (517) 

which can also be obtained by the straightforward substitution c = 2C into Eq. (I5.16p . 



6 Geometrical formulation 

Special relativistic kinematics can be given a geometric interpretation introducing the Minkowski 
quadratic form 

rj(dt, dx) = dt 2 - dx 2 /c 2 , (6.1) 

which is invariant under Lorentz transformations of the coordinates. We now ask what is the 
most general function of dt and dx that is invariant under the transformations found in Sec. [2j 
That is, we look for a function ( such that 

C(di, dx) = C(dt, da;) , \/v e J . (6.2) 

Differentiating Eq. (16.21) with respect to v and evaluating the result for v = we find 

This is a particular case of the partial differential equation solved in Appendix [UJ corresponding 
to: X = dt; Y = dx; A = 1/k - 6/2; B = a; C = 1; D = 1/k + 6/2; F = (; V = tp\ Q = h. 
Hence, 

C(dt, dx) = e - 2 M<Wdt)/ K ^2 + hdtdx _ a dx 2) ^ (g ^ 

where the expression (12.271) for (p has been used, is the extension of the quadratic form (16.11) . 
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Anisotropic kinematics can then be interpreted in terms of a spacetime geometry by intro- 
ducing the pseudo-norm 

m = e -*W)/« ((£0)2 + H o e _ fl ( ^i )2) i/2 j (6 5) 

where £ is a generic vector in spacetime. The most general spacetime structure compatible with 
the relativistic framework considered in this paper is thus a pseudo-Finslerian, not a pseudo- 
Riemannian F^l one. The pseudo-Finslerian character, associated with the exponential factor 
in Eqs. (16. 4p and (16. 5p . is unavoidable when one wants to allow for a mechanical anisotropy, 
through the parameter k. It is interesting to note that the other kind of anisotropy, linked 
to the parameter e, does not alter the pseudo-Riemannian character of spacetime. This can 
be exemplified in the physically interesting case of anisotropic Einstein's kinematics, where 
Eq. (16 .4p become^ 

C(dt, dx) = (dt + dxlc + f- clK (dt - dx/c_) 1+c/K . (6.6) 

When \k\ = +oo but s ^ 0, this is just the Minkowski quadratic form (16. ip in non-Lorentzian 
rectilinear coordinates, defined operationally through a non-standard synchronisation proce- 
dure. The fact that the value of e does not alter the spacetime structure agrees with the thesis 
that the anisotropy associated with such a parameter is an artifact of a convention rather than 
a physical feature (see the discussion in Sec. 12.2. 2p . Indeed, one can always eliminate such 
an anisotropy by synchronising clocks according to the Einstein procedure, which amounts to 
choosing Lorentzian coordinates in spacetime. On the other hand, a finite value of k implies 
true physical effects (anisotropic time dilation and length contraction, modified dispersion rela- 
tion, ...), that cannot be gauged away by a stipulation, just as the pseudo-Finslerian character 
of spacetime cannot be reduced to a pseudo-Riemannian one by a suitable choice of coordinates. 

For a particle tracing out a differentiable worldline x a (9) in spacetimel^l where 9 is an 
arbitrary parameter, we can define the tangent vector with components U a = dx a /d8, and the 
proper time r such that dr = F(U) d6. Of course, U 1 = uU°, so by a comparison between 
Eqs. (14.11) and (16.51) we note that for \k\ ^ c, 

L <"> di = -TJ^W F(U) d9+ (( P "- m T0^W) 1,1 ~ ( £ » " W^W) U °) M ■ 

n (6 ' 7) 

Hence, apart from a constant additive termJ^J the action turns out to be proportional to the 
particle proper time, as it happens in ordinary (isotropic) special relativity. 

Two comments are in order, of opposite flavour. First, we see that the Newtonian La- 
grangian mu 2 /2, in spite of its appearance, is related to spacetime geometry. For it is the 
limit, when \k\ —>■ +oo, of the Lagrangian (I4.12p . which is the particular case of Eq. (16. 7p 



24 The difference between Finslerian [29] and pseudo-Finslerian [30 structures on a manifold is the same as 
between Riemannian and pseudo-Riemannian ones. 

25 In the case c+ = c_ (i.e., e = 0), this expression was also considered in Refs. [2"Tl |2"21 |3"T] . 
26 The indices a, b, ... run from to 1. 

27 Which vanishes when the rest momentum and energy are given by the expressions (|4.7I) and (|4.8[) . 
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corresponding to the norm 

^(0 = exp (6.8) 

and to c = +00. Second, although it is customary in the literature about relativity to choose a 
Lagrangian such that the corresponding action turns out to be related to the proper time, this 
should not be taken as a dogma. In the exceptional cases with k = ±c, the pseudo-norm is 



e-[-- £ -)e 

K C 



(6.9) 



so F{U) d8 = (1 — (1/k — s/c) u) dt, which does not contribute to the equation of motion. 
Hence, not only Eq. (\6.7\i does not apply to these cases, but no term proportional to proper 
time in the action leads to a non-trivial dynamics — in fact, Eq. (I4.15P does not appear to have 
a geometrical interpretation. 

Equation ( 16. 71) allows one to define a Lagrangian in spacetime, 

^ = -TJ^W F[U) + (( P ° " m W^W) Ul - i E ° - W^w) U ' 

(6.10) 

so one can define a momentum one- form in spacetime, with components p a = dC/dU a . It is 
easy to see that p = —E(u) and p x = p(u). Note that, since £ is a homogeneous function of 
degree 1, the components p a are independent of the parametrisation of the particle worldline. 
It is also worth stressing that the relation between p a and U a is not trivial. Indeed, on defining 
the quantities po = —E and p\ = p as we did in Sec. El one finds 

dF(U) _ m g ab (U)U b , , 

Pa 1/c 2 -1/k 2 8U a 1/c 2 -1/k 2 F(U) ' K ' ] 

where in the last equality we have used Euler's theorem for the function dF(U) 2 /dU a , which 
is homogeneous of degree 1, and the definition [301 E2] 



It is now possible, introducing coefficients g? b {p) such thato g? b (p)gb c (U) = 5 a c , to write the 
dispersion relations of Sec. 15.11 in the covariant form [32J 



36/^ ^ m 



® ? - ft= (I/,? r !/«»)» ■ (6,13) 



where the relation g a b(U) U a U b = F(U) 2 , which follows from the fact that F 2 is homogeneous 
of degree 2, has been used. 



28 The g ah are homogeneous functions of degree of their arguments (as the g ab are), so their definition is 
insensitive to the coefficients in Eq. (|6.11|) . 
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7 Comments 

In this paper we have developed the foundations of particle dynamics for theories that are 
compatible with the relativity principle and include a mechanical anisotropy. We have seen 
that anisotropy can be taken into account through two parameters, one of which (e) has a 
merely conventional status, while the other (k) accounts for real physical effects. It may then 
seem odd, that the parameter e enters in the expressions of physically relevant quantities such 
as momentum and energy; for this appears to offer an opportunity to test conventional isotropy 
by means of dynamical experiments, which should of course be impossible. In fact, the value of 
e cannot be dynamically determined, independently of kinematical considerations, because any 
such experiment requires that one could also determine velocities, which again presupposes a 
synchronisation procedure. Thus, the presence of e into equations such as H4.5[) and (14. 6[) is no 
longer surprising, its role being analogous to that played by a conversion factor. 

On the contrary, the value of n can be determined experimentally. For example, one might 
study momentum conservation in a collision. Having synchronised clocks with some procedure 
(corresponding to a value for e) and defined velocities correspondingly, one can look for the 
value of n that provides a best fit when used in Eq. (14.51) . Another possibility, less direct but 
perhaps more practicable, is to place constraints on the deviations from standard dispersion 
relations, as it is now fashionable to do within the context of the so-called "quantum gravity 
phenomenology" [55] . 

We have deliberately restricted our investigation to those cases in which time and space are 
homogeneous, and space is geometrically homogeneous and isotropic; the only manifestations 
of anisotropy being mechanical. This choice is motivated by the desire to set up a kinematical 
framework based on reference frames of the type usually considered in Newtonian mechan- 
ics and special relativity. Of course, it would be interesting to generalise mechanics further, 
perhaps even contemplating situations in which space and time are inhomogeneous, possibly 
at short scales. (This is what several people believe might happen in a more sophisticated 
theory of spacetime structure.) However, such generalisations must necessarily proceed along 
lines that differ from those of the present paper. For example, one might explore anisotropic 
deformations of de Sitter relativity [51] , or try to include anisotropy at the level of the algebra 
of generators [55] . 

Recently, motivated by experimental data about ultra- high-energy cosmic rays [55], there 
have been some suggestions to modify the standard dispersion relations in such a way that 
the relativity principle is preserved, allowing at the same time for the existence of an invariant 
energy scale. Such modifications have been developed essentially within two theoretical frame- 
works: Nonlinear representations of the Lorentz group [57], and the so-called K-deformations 
of it [55] (see also [39]). In both cases, the emphasis is on energy- momentum space, while the 
role of configuration variables, such as position and velocity, remains somewhat unclear [40] . 
Of course, the corresponding dispersion relations differ from those found in Sec. [5j hence, these 
frameworks are incompatible with at least one of the hypotheses on time and space made in 
the present derivation. (It is also possible that they are not fundamental as suggested re- 
cently |41j.) Another, less fancy, origin for modified dispersion relations could simply be a 
breaking of Lorentz invariance, perhaps combined with small-scale inhomogeneity. This hap- 
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pens, e.g., if one replaces space by a discrete structure, in which case Lorentz invariance is 
recovered only at low momenta, which do not probe the lattice structure of the background 
and thus the fundamental asymmetry between space and time. A similar behaviour is exhibited 
by condensed matter models, such as Bose-Einstein condensates [42J. 

Using the equivalence principle in order to argue that gravitational effects can be locally 
gauged away, one can construct a theory of curved pseudo-Finslerian spacetime that locally 
reduces to the one of Eq. f 1 6 . 6 j) . thus obtaining an extension of general relativity. One might 
even construct a theory in which the amount of anisotropy changes from place to place, e.g., 
replacing the constant parameter k by a field [22J. In this case, it is possible to envisage 
situations in which spacetime here and now is almost Lorentzian, whereas elsewhere (perhaps 
in regions of strong gravity, or in the early universe) it is highly non-pseudo-Riemannian. 
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Appendices 

A Alternative derivation of the transformation law 

Replacing Eq. (12. 3ft and the corresponding expressions for A(u) and A($(u,t>)) into Eq. (12.7ft 
we find, after elementary manipulations: 

A(v)A(u) (1 +u£(v)) = A{<f>{u, v)) ; (A.l) 

V + UT](v) , . . 

= ird55 ; (A - 2) 

+ ((v)ri(u) = (l + u£(v))Z(<S>(u,v)) ; (A3) 

v£(u) +rj(u)n(v) = (l + u£(v))ri($(u,v)) . (A.4) 

From Eq. (1A.2j) we immediately get the expression (12.271) for <p(u). Differentiating Eqs. (1A.3j) 
and (IA.4I) with respect to v and evaluating the result for v = 0, one finds: 

arj(u) = au£(u) + (l + bu — au 2 ) £'(«) ; (A. 5) 

£(w) + bt]{u) = auj]{u) + (l + bu — au 2 ) rj'(u) . (A. 6) 

Replacing £(«) from Eq. (IA.6j) into Eq. (IA.5I) . we obtain rj"(u) = 0, which is trivially integrated 
with the conditions 7/(0) = 1, t/(0) = b, to obtain Eq. (I2.22p . Equation flA.6[) gives then 
immediately the expression (12.211) for £(«). 
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The function A is still undetermined, but we can now find its general form by replacing the 
expression (12.21 j) into Eq. (lA.lj) . then differentiating the resulting equation with respect to v, 



and finally evaluating the result for v = 0. One finds, remembering that A(0) = 1, 

A'(0)A(u) + au A{u) = A'{u)(f{u) . (A.7) 
This equation can be immediately integrated to obtain Eq. (12.291) . 

B On the mass-energy relation 

Consider, in a reference frame /C, two particles with masses mi and m 2 , and velocities U\ and 
u 2 , such that 

p(mi, ui) + p{m 2 , u 2 ) = p{M, 0) = p {M) (B.l) 

for some M, function of mi, m 2 , u\, u 2 . Thus, K, can be regarded as the centre-of-momentum 
frame for the system!^! and M as its total mass. In a reference frame /C, in which K, moves 
with velocity v, we have 

p(mi, $(«i, v)) + p(m 2 , $(w 2 , v)) = p(M, v) . (B.2) 



The arbitrariness of v allows us to find the value of M. Indeed, differentiating Eq. (1B.2j) with 
respect to v, and setting v = 0, we get 

p'{ mi , + p'(m 2 , u 2 )cp(u 2 ) = p'{M, 0) = L"(M, 0) = M , (B.3) 

where a prime denotes the derivative with respect to velocity. Remembering now Eq. (13.351) . 
Eq. <KW\ with A = 1, and Eqs. (13T29D and (EQ2l . we find 

p'(m, tt)</?(w) = — (p(m, «) — po(m) — // (E(m, u) — E (m))) . (B.4) 

Note that, whereas in Einstein's dynamics (p = Eu/c 2 , po — 0, fi — 0) the right-hand side of 
Eq. (IB.4[) turns out to be proportional to E(m, u), this is not the case in general. 



Let us restrict ourselves, from now on, to the anisotropic dynamics of Sec. 14.11 It is conve- 
nient to define the function 

I» := e~ h ( u)/K <p(u)- 1/2 . (B.5) 
Then, the momentum and energy in Eqs. (14. 5p and (14. 6 p can be rewritten as: 

m ( 1 — e 2 \ \xm 

P(m, u) = 1/c2 _ 1/k2 ^ + — uj T(u) + poH - 1/c2 _ 1/k2 ; (B.6) 

TYl Til 

E(m, u) = 1/c2 _ 1/k2 (! + (// + 2e/c) u) T(u) + E (m) - — _ — . (B.7) 



29 Whcn e = and \k\ = +oo, the centre-of-momentum frame can equivalently be defined as the one where 
total momentum vanishes. This is not the appropriate characterisation when anisotropy is present in some 
form, as one can realise considering a situation in which the system is made of a single particle. Also, note that 
defining a centre-of-mass frame is problematic in mechanics where there is no absolute time |43j . 
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It is now easy, remembering Eq. (13.301) . to see that p'(m,u)<p(u) = mT(u), so 

M = mi r(ui) + m 2 Y{u 2 ) . (B.8) 

Assume now that Eq. (14.71) holds. It is then easy to see that Eq. (IB. II) amounts to 

mi r(«i) iti + m 2 T(u 2 ) u 2 = . (B.9) 

Using this result in computing the total energy of the system in the centre-of-momentum 
reference frame JC, and assuming also the validity of Eq. (14. 8p . one finds 

M 

E(m u Mi) + E{m 2 , u 2 ) = - 2 . (B.10) 

Hence, Eq. (14.81) holds also for the system as a whole, not only for the individual particles. Of 
course, this conclusion can be straightforwardly generalised to an arbitrary number of particles. 

Note that this argument goes through only if po and Eq are of the form given by Eqs. (14. 7\i 
and (14.81) . In fact, although the result concerns the notion of rest energy, its validity requires 
also that there be a rest momentum p = pE Q . This is not an option: If one accepts the 
existence of a rest energy, such a rest momentum must also exist, for consistency, as can be 
seen by the following argument. In the centre-of-mass reference frame /C, 

E(m 1 ,u 1 )+E(m 2 ,u 2 ) = E (M) . (B.ll) 

In the arbitrary frame /C, 

E(m x , $(mi, v)) + £(7712, $(«2, v)) = E(M, v) . (B.12) 

Differentiating Eq. (IB.12P and setting v = 0, we find 

E'{ mi , mi)</?(mi) + E'{m 2 , m 2 Mm 2 ) = E'{M, 0) = , (B.13) 

where the right-hand side can be evaluated, e.g., using Eq. (13.351) and the fact that p'(M, 0) = 
M. Using now Eq. (13~TS|) with A = 1, Eqs. <Km and (I5^2]l . and finally Eqs. ([B3]) and (IB~lTD . 
Eq. (IB. 131) becomes 



p (M) -nE (M) -po(mi) - p (m 2 ) + p, E ( mi ) + u, E {m 2 ) = . (B.14) 

Remembering now that both po and Eq must be proportional to the mass (see end of Sec. 13.41) . 
setting po = a Eq, and using Eq. (1B.8|) . we find easily that a = p. 



C Solving a partial differential equation 

Consider the first order linear partial differential equation 



(AX + BY)^ + (CX + DY)^ = 0, (C.l) 
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where A, B, C, D are constantsjfj Defining the new variable Z := Y/X and the function 
G(X, Z) := F(X,XZ), Eq. (JSHJ can be rewritten as 



1 dG *(Z) dG_ 

2 dX + -V(Z) + {A + D)dZ~ ' 1 J 

where 

:= -fiZ 2 - (A- J D)Z + C . (C.3) 

Defining 

6(2) := / (C4) 



and a new variable through 



dW := ^ - (A + D)dQ , (C.5) 



Eq. (1C.2|) can be immediately rewritten in the form 



l Y dG dG n (r^ 
2 X dX~dW = °' (C - 6) 

where G now denotes, with a little abuse of notation, the function G when Z is expressed 
in terms of W. It is evident that G can be any arbitrary function of lnX 2 + W, so F is an 
arbitrary function of 

X 2 ^{Y/X)e^ A+D ^ Y ' x ^ . (C.7) 
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